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Abstract 

We consider random primordial magnetic fields and discuss their dissipation, 
coherence length Lq, scaling behaviour and constraints implied by the primoridal 
nucleosynthesis. Such magnetic fields could excite the right-helicity states of 
Dirac neutrinos, with adverse consequences for nucleosynthesis. We present 
solutions to the spin kinetic equation of a Dirac neutrino traversing a random 
magnetic field in the cases of large and small Lq, taking also into account elastic 
collisions. Depending on the scaling behaviour and on the magnetic coherence 
length, the lower limit on the neutrino magnetic moment thus obtained could 
be as severe as lO^^^fj^B- 
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1 Introduction 



The phase transitions of the very early universe may have generated primordial mag- 
netic fields, which could play an important role in cosmology. Large primordial fields 
can survive to this day and provide a seed field for the galactic dynamo mechanism 
[0] which amplifies the seed to produce the observed galactic magnetic fields. Several 
suggestions have been made as to the possible mechanisms which could produce large 
primordial fields f^. For instance, large fiuctuations at the electroweak phase transi- 
tion might be responsible 0, and it has been argued that after a proper statistical 
averaging one could actually obtain seed fields of the required magnitude about 
10~^^ G. A more exotic possibility which also seems to work is based on the obser- 
vation that in Yang-Mills theories the vacuum may have a permanent magnetic field, 
which is imprinted on the comoving plasma already at the GUT scale [§. 

Primordial nucleosynthesis is sensitive to magnetic fields, which modify both the 
Hubble rate and the rates of the reactions that help to build up the light elements 0]. 
Primordial nucleosynthesis may also be affected in another way p| provided neutrinos 
are Dirac particles. In that case the right-helicity component of the neutrino can be 
excited and brought into thermal equilibrium by scattering of the left-helicity neutrinos 
off the magnetic field, thus changing the effective number of degrees of freedom. 

In a recent estimate of the primordial helium abundance the allowed number of 
the additional neutrino species was reduced to AA''^ = A^^i, — 3 ~ 0.1. There is some 
uncertainty in this estimate due to the unkown magnitude of systematic errors in the 
observed abundances, but even allowing for very conservative systematic errors, AA^j, 
is definitely less than 1. Thus the coupling of a Dirac neutrino to a primordial magnetic 
field should be weak enough not to equilibrate the right-helicity states below the QCD 
phase transition temperatures. This gives rise to a bound on the combination /i^-B, 
where /i,^ is the magnetic moment of the neutrino. 

The right-handed neutrino production rate is proportional to the neutrino helicity 
flip probability, which may be calculated by considering neutrino spin rotation in 
a medium with an external magnetic field. The complicated time evolution of the 
neutrino spin can be described in terms of a relativistic kinetic equation (RKE), which 
has been derived in and was extended to account for the elastic collisions in 
The helicity rotation of a light Dirac neutrino traversing a magnetic field is determined 
by forward scattering off the field. Simultaneously, the neutrino interacts with all the 
particles in the plasma via reactions that for light neutrinos can be taken to conserve 
helicity. An essential feature in deriving the RKE is averaging over the (random) 
magnetic field, the procedure of which we improve in the present paper in order to 



1 



discuss also small scale fields. In what follows, we shall limit ourselves to the Standard 
Model but assume that the neutrinos have Dirac masses. 

If we neglect the magnetic moment of the Standard Model Dirac neutrino, in the 
ultrarelativistic limit the dispersion relation, in hot plasma and in the presence of a 
magnetic field, reads E{q) ^ q + V. Here V is the neutrino interaction potential which 
is determined by the neutrino forward scattering amplitude and which consists of two 
parts [0: 

V = V^'"'^'^^ + v^''^^"-^\ (1) 

yiaxtai) j^cludes coutributious which are due to the fact that in the presence of a 
magnetic field the charged background is actually magnetized. This turns out to be 
an important effect. The interaction potential V^'"^'^^ is generated by the mean vector 



current ~ {'ipa7fj.'4'a)o, and at the temperature rUe <^ T <^ Tqcd it reads [jT3 



yivec) ^ 3 4 ^ 10^2° (^^)^ MeV . (2) 

yiaxtai) generated by the mean axial current of the magnetized leptons ~ {'4^n^lb'^i)o 
and is given by |]l2l 

V<°"°" = f.e„^ + ^(B=-^), (3) 

where the quantity /Xg// is defined byQ 

eG'F(-2cA)Tln2 ,0, , / T n 

Mf = ^'J 6 X 10--(-2c,);.,(^) , (4) 

and ca = =f0.5 is the axial constant in the weak lepton current (upper sign for Ue, 
lower for z/^,t)- Here fiB = e/2me is the Bohr magneton. In what follows we may 
safely neglect the last term in Eq. (||). 

The size of the random magnetic field domain Lq influences crucially the neutrino 
spin behaviour. For large-scale random magnetic field with Lq ^ V^~^ the dominant 
mode of evolution is spin oscillation. For small scale magnetic fields (which is perhaps 
a more realistic alternative in the early universe) there also appears aperiodic spin 
motion which effectively converts to uji. We discuss Lq, dissipation, conductivity 
and general constraints on primordial magnetic fields in section 2. In section 3 we 
derive a cosmological limit on the neutrino magnetic moment in the case of a large 
scale magnetic field, taking into account elastic collisions and improving on previous 
treatments. Section 4 introduces an averaging procedure which is suitable for small 



^Note that this effective magnetic moment has no relation with the anomalous neutrino magnetic 
moment. 
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scale magnetic fields, and we provide a limit on the neutrino magnetic moment also 
in this case. Section 5 contains our comments and a discussion of the meaning of the 
results. 



2 Direct constraints on magnetic fields 

The early universe is an excellent conductor, and a primordial magnetic field, if such 
existed at any time, is imprinted on the comoving plasma which will retain the field. 
The magnetic flux is conserved so that the magnetic field scales with the expansion 
of the universe as i? ~ R~^. The strength of a random magnetic field depends also 
on how it scales over physical distances. If the field were to perform random walk in 
3-d volume, the scaling would be i? ~ N''^^"^, where N is the number of steps. An 
argument based on the statistical independence of conserved flux elements gives rise 
to a 5 ~ A^~^ scaling. Whether the magnetic flux actually is completely uncorrelated 
in the neighbouring unit cells is an unsolved issue. In |Q it was argued that if the 
magnetic field domains are uncorrelated, a proper statistical averaging produces a 
field with {B) = and a root-mean-square field Brms = ~ A^~^/^. In the present 
paper we take a more phenomenological view and assume merely that there exists a 
random magnetic field with the scaling 

Rn \ f Lq\^ 



where L and Lq are two comoving physical scales, and p is essentially an unknown 
parameter; we shall focus mainly on the choices p = 1/2, 1, 3/2. We shall view 
B{L,t) in Eq. (|^) as a root-mean-squared field with (Bj) = 0. We shall also assume 
isotropy so that {BiBj) = for i 7^ j. If we further assume that the magnetic domains 
are uncorrelated we may write 

(A(x)B,(y)) = (2A)-i5,/=')(x-y) , (6) 

where the length A is determined by the domain size Lq and the rms field (|^) at the 
horizon scale L = Ih{T): 



= 7r(3-2p) ^'-^^ = ^^^^°' P^3/2, 

= |(ln^)i?L.(^ = M^^ P = 3/2. (7) 
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We may prove Eq. (|^) by recalling the general fluctuation theory formula^ [|I4| 

(5,(k)5,(k')) = {2nf{B\6^''\k + k'){6i, - hk) ■ (8) 

If we substitute the Fourier transform -Bi(x) = / (i'^A; exp(zkx)i?j(k)/(27r)'^ into the 
Anzatz (H) and use the rule Eq. (H), we obtain a constant spectrum 

(B')k = \^ = const , (9) 

where the factor 3/2 arises because of a different delta-function implementation in Eq. 
(H) as compared with Eq. . 

Now we can express the parameter A in Eq. (^) via the field in Eq. (j^). Let us 
first note that the mean magnetic field energy within a horizon volume, 



depends^ on a cut-off wave number kmax'- 

Vh ^Tq^ Jlo ^ r ^ IbAvr^ 

Physically k^ax = Stt/Lq corresponds to a minimum size of the inhomogeneity, and is 
given as in [10| by 

3, _ Jk'd'k{B'h/{2nr 

r max /T39\ ■ V / 

Here the factor 3/5 arises by using Eqs. @ and ([l0|). 

Substituting Eq. (|T2p into Eq. (|TTp we arrive at the relations Eq. (|^ between A, 
the domain size Lq and the field Brms at the horizon scale Ih'- 

and similarly for p = 3/2. 

Let us now turn to the constraints on the domain size Lq and on the magnetic field 
strength in the plasma of the early universe. We begin by discussing the dissipation 
of primordial magnetic fields. 

^For real fields we use the + sign in the (5-function argument instead of the usual - sign for the 
conjugated magnetic field components. 

^The second equality here has the sense of the mean field energy density i(B^)x=o = 
1 1 (27r)^ (^^)k for fluctuating fields jl5|. The first equality in Eq. (|l^) has this meaning as a 
definition. Of course, for any index p in Eq. (|) the magnetic field does not affect the expansion of 
the universe because pB P-y ^ T'^- We use here the Gauss units in the which for the uniform field 
(B^)k = i?Q(27r)'^(5(^)(k) the magnetic field energy density is given by ps = 
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Conductivity in the early universe, although large, is not infinite, and accordingly 
there will be some ohmic dissipation, starting from small length scales, Lq < L^iss . 
The dissipation scale is given by 



-'diss 



(14) 



where t is the age of the universe and cxc is conductivity, which for isotropic relativistic 
electron gas reads 

o;^ T 

<yc = -r^ = • (15) 

AiraUp iraLr, 

2 T /T^2 



Here Up = y Anarie/T is the plasma frequency and a ~ na Lc/T is the collision cross 
section and Lc ~ 3 is the Coulomb logarithm. We then find that the dissipation length 
below Tqcd is given by 

rrel n 1 -l/4/MeVx3/2 

L^ss = 0.1(7, ^ cm . (16) 



The finite isotropic conductivity (|15D refiects the properties of the plasma at large 
scales when we may neglect the random magnetic field infiuence upon relativistic 
plasma. We may thus use the estimate (0) to obtain a conservative lower bound on 
the domain size Lq. 

It is worth noting that within a domain volume with a uniform magnetic field where 
plasma is highly anisotropic, the dissipation length is greatly reduced with respect to 
Eq. (17) because of substantial increase of conductivity, i.e. a strong magnetic field 
tends to sustain itself within a given domain. Indeed, for relativistic electron-positron 
plasma in a strong uniform magnetic field, which is much larger than the quantizing 
magnetic field Be = m^/e = 4.41 x 10^^ G, conductivity is given by 



Here a ~ 4:C'^TTa'^{Bc/B)/T'^ is the e~^e collision cross section [16 for the electrons 



and the positrons that occupy the lowest Landau level, and Ce ~ 0.577 is the Euler 
constant. The dissipation length within a domain with a strong uniform field is given 

by 

_ 3.5x10-^ cm ^MeV^3/2 .5,. 1/2 

Ld^ss = -T7i [-^) X i^j • (18) 

y* 

Let us note that within one domain with a uniform magnetic field the ratio B/Bc can 
be interchanged with the temperature ratio, B / B^ ^ (T/MeV)^. The dependence of 
the magnetic field on the temperature refiects the magnetic fiux conservation. Here we 
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neglected the small contribution to the plasma conductivity produced by the collision- 
less one- photon processes in a strong external magnetic field (7 e~^e~, 6^7). 
This conductivity is calculated from the imaginary part of the polarization tensor in 
magnetic field [jl^ and can be estimated as 



^one-photon ^ C^(^Agc)m,^ ^ ^ ^^g^ 

where cTc is given by Eq. (p!7|). 

Comparing the relativistic dissipation length ([161) with the analogous result for 
anisotropic plasma ([T8|) we find that one can neglect the magnetic field dissipation at 
small scales since domains with sizes Lq{T) > LJJ^^^(T) ^ L'^^l survive. 

In the non-relativistic case, which is relevant for the recombination time t ~ 10^^ sec 
{Tree ~ 0.4 eV) and for the BBN time t ~ 1 min (Tbbn ~ 0.1 MeV), one obtains a 
different expression for the dissipation length, 

ir,-"' = 0.l9.-"'(^)"*cm, (20) 



where we have used the conductivity of non-relativistic isotropic hydrogen plasma p8 



= —T^^ 172 ■ (21) 



The difference between this and the relativistic expression (|T^) follows from: (i) the 
Debye screening for the forward ee or ep scattering in isotropic non-relativistic plasma 
with the Coulomb logarithm Lc ~ 10, and (ii) the fact that in the non-relativistic case 



the velocity of the electrons is the thermal velocity vt^ = ^jT/rrii. 7^ 1. 

The dissipation length (^) at the recombination time is about ~ 10^*^ cm which 
translates to the BBN time t ~ 1 min {Tbbn ~ 0.1 MeV) as L^'''{Tbbn) = 10^° cm x 
{Trec/TBBN) — 4 X 10^ cm. The domains which at that time are larger than this 
scale, Lq>L^^"-{Tbbn), survive after the recombination time, and such a relic field 
could be a seed field for the dynamo enhancement of the galactic magnetic field. This 
requirement can be rewritten as (see also 0) 

Lo > L™" = 10^ cm X (^) . (22) 

Let us note that the local dissipation length calculated at the same BBN time from 
Eq. ( ppp is sig nificantly smaller, L'^Z'^'^Tbbn) 

~ 5 cm, i.e. domains with sizes within 
the region 5 cv!i<Lq{Tbbn)'^L'^^'^ {Tbbn) ~ 10^ cm continue to expand after the BBN 
time but they dissipate before the recombination time and do not contribute to the 
relic seeding of galactic magnetic field. 
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For instance, the relativistic dissipation length ( [T6| ) at the temperature correspond- 
ing to the electron neutrino decoupling with matter, T ~ ~ 2 MeV, is even smaller. 



L^t(T ~ 2 MeV)~ 10-2 cm, (23) 

or ~ 0.3 cm after its expansion till the BBN time. Therefore, all the domains 
that have survived in the hot plasma but have at Tbbn a size within the range 
0.3 cm<Lo(TeBAr)<5 cm, in fact dissipate before Tbbn- On the other hand, all the 
domains that survive after the nucleosynthesis but have a size less than 10^ cm at 
Tbbn dissipate before the recombination time. 

Accounting for the expansion of the domains, we find that if 

Lo(T)>L— ' = 10-2 cm x (^) , (24) 

such domains are guaranteed to survive at any relativistic temperature; however, they 
dissipate before the BBN time if their scale Lq{Tbbn) was less than L^°^~'^'^\Tbbn) = 
5 cm. 

Primordial nucleosynthesis considerations can further be used to constrain the 
strength of the magnetic field. If large enough, a magnetic field would affect the 
expansion rate of the universe as well as the rates of the various reactions that are 
involved in building up the abundances of light elements. Including these effects in a 
nucleosynthesis code Cheng et al. find in a recent study an upper bound 

5ns = 5(10^ cm, 1 min) < 10"G. (25) 

This is a very useful constraint, as will become evident shortly. 

The energy density in the magnetic field within the causal horizon is given by 

where Vh is the causal horizon volume. (At Tew — 100 GeV the causal horizon size 
is Ih — 1-4 cm). A reasonable requirement is that the energy of the magnetic field 
fluctuations should not exceed the free electromagnetic energy, so that pB<T^. This 
is not in disagreement with the nucleosynthesis constraint Eq. (^) at any scale below 
Tew- If we assume, as is natural, that the random field fiuctuations are bounded by 
B{LQ,t)<0{T'^), then pb'^T^ for Iq < lu for all p under consideration. 

If we assume that the observed galactic magnetic fields are due to a large scale 
dynamo effect m, seeded by the primordial field, then we may set a lower bound on 
the strength of the cosmological seed field -Bseed (this is further enhanced by a factor of 
10'* by the collapse of the protogalaxy) [|1^. Numerical simulations of the field growth 
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appear to imply that the seed field must be sufficiently large for the dynamo to work 
0, and one should require that at t ~ 10^°yr and at the intergalactic scale L^ai = 100 
kpc, -Bseed = -B(-LGai, lO^^yr) > lO^^^G. This bound can be made even tighter provided 
we assume that the observed field reversal between the Orion and Sagittarius arms 
20[] is related to the dynamo. To produce a field reversal, the seed field should have 



been relatively strong |21|, i?soed^lO^^^G. Such reversal has only been observed in the 
Milky Way and may thus not be a generic feature. 
From Eq. (p5|) it follows that 



B... = Bns (^) i^X" (l^HH^y < 1.8 X 2.4' X 10-'-"'G. (27) 

V t* / \tnow/ \ -t^Gal / 

where we have assumed that the change to matter dominated universe takes place 
at ~ 8750 yr, corresponding to Qh"^ ~ 0.4. Therefore, the case p = 3/2, which 
translates into the bound -Bsoed < 2 x lO^^^G, cannot provide a seed field large enough 
for the galactic dynamo, no matter what mechanism generated the field in the first 
place. If p = 1 we obtain -Bseed < 4 x 10~^^ G but given the large theoretical uncer- 
tainties, this case could still be compatible with the galactic dynamo. Finally, if the 
observed field reversals in the Milky Way are really related to the dynamo, then only 
the p = 1/2 case remains marginally compatible. 



It has been argued ||22| that in magnetic fields of the order of -B ~ 10^^ G the 
neutron becomes stable against /5-decay and that for somewhat larger fields proton 
becomes unstable to a decay into a neutron. If B<T^ at all scales, this effect would 
not be important for nucleosynthesis. It is also irrelevant for the problem of excitation 
of the wrong-helicity neutrino states because neutrino spin ffip is determined by the 
field strength at Tqcd and at scales of the order of the weak collision length L^. 
Indeed, we may write the the nucleosynthesis bound (^Sj) as 



B{Lw, Tqcd)<8 X 10^^ X 13-PG < 2 x lO^^G (28) 

Here we have used the estimate ^ 1.6 x 10~^ m. Considering the inherent uncer- 
tainties, we may then well adopt the latter figure in Eq. (^Sj) as a conservative bound 
at Tqcd for all p. 

Finally, note that the evolution of a magnetic field in an expanding Universe is 
determined by MHD which should explicitly yield the domain size Lq as well as the 
topological index p in Eq. (^, both of which are crucial for the neutrino propagation 
in a medium with a random field. We shall use the phenomenological model (^ for 
case of uncorrelated 3-d fields in Eq. (P) with an arbitrary index p. 
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3 Neutrino propagation in medium with large-scale 
random magnetic field 

Neutrino spin-flip in a magnetic field can affect the Big Bang nucleosynthesis of light 
elements because of the appearance of an additional gravitating relativistic component 
in the plasma. The nucleosynthesis limit on the extra degrees of freedom at the time of 
nucleosynthesis is 0, in units of relativistic two-component neutrinos, AA''j,<0.1 (see 
also l2^). Indeed, the wrong-helicity neutrinos will be abundant at the time when the 



neutron-to-proton ratio freezes at T ~ 0.7 MeV, violating the nucleosynthesis bound, 
unless they decouple before the QCD phase transition. Then their relative number 
densities will be diluted to acceptable levels by the heating of the particles still in 
equilibrium. 

In the early universe the production rate of the wrong-helicity neutrinos is given 
by Tl^r = {Pi^j^^^j^)T]^ , where T]^ is the total weak collision rate, and in the absence 
of collisions the averaged conversion probability reads 

(P^.^^.) = (29) 

where H± = 2fii,B± is the field perpendicular to the neutrino propagation, fii, is the 
neutrino magnetic moment and 

uj^ = 2Hl + V^ + QL^y^ (30) 

is the spin rotation frequency. Note that it depends both on the value and scale of 
the random field. In a uniform constant field one would find uj^ = H\ + V"^ in the 
expression (|29|). In a realistic case V S> H±. 

To find out Tl^r, let us consider the electron neutrino for definiteness. Then at 
T = Tqcd one should include the following processes: z/cZ/c , i^ai^a, T-'el^ — ^ 

z/^/i"*"; z/eZ/^ — fi^e^. Here the notation is: a = e, /i, r; / = e, /i ; q = u, d. Adding 
up all these processes we find the thermally averaged electron neutrino collision rates 
to be 

^v/{Tqic-d) ~ SOGITqcd; r^(TQCD) ~ 2G|Tqcj3. (31) 

Here we have neglected the decays and inverse decays, whose contributions are small. 
The main contribution to Eq. (^) comes from the charged current processes involving 
quarks. For and Vr the rate is slightly different. 

Requiring Tl^r<H at T = Tqcd where H = ^Sii^g^T) /90T^ /Mpi is the Hubble 
parameter, one obtains a constraint on the product of the Dirac neutrino magnetic 
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moment and the mean squared random field B = (5^)^/^ [§. Here we adopt the value 
5'*(^qcd) — 63, which includes also the effects due to the non-relativistic species [|4 . 
The limit thus obtained does not depend on any model for the random field, and one 
finds that 

^^MTQCB,Lw)<4 X 10Vi.(^^^)'/' G . (32) 

The result (^) is true in collisionless regime if the random domain size Lq is larger 
than the neutrino spin oscillation length lose — with 

lose = IQ-^lHiT/MeV)-^ , (33) 

where In ~ Mpi/T'^ is the horizon length. For small domains, Lq <^ lose, neutrino spin 
cannot follow the random direction of the magnetic field and the neutrino spin rotation 
effectively ceases The limit ( |32D is slightly different from the one presented in 



because of the inclusion of the quarks in the total collision rate. 

If we take into account elastic neutrino collisions with charged particles in the back- 
ground plasma, the analogous consideration starts from the RKE for the 2;-component 
of the neutrino spin Sz{t) = 2P^^^^^(t) — 1. This is an integro-differential equation 



and given by fTT 
dt 



2 /o* exp[- X u^{h)dt^] X [H_{t)H+{t,) exp(z //^ V{t,)dh) 

+H4t)H4h) exp{-tJl V{h)dh)]S,{h)dh , (34) 



where H±{t) = jiu{Bx{t) ± iBy{t)) depends on the transversal field components only. 
The collision frequency z^±(t) is approximately given by the weak elastic collision rate 
Fp^. In the absence of inelastic collisions the relaxation of the transversal spin com- 
ponents is completely determined by Eq. (0). If V and v±_ are slowly varying, one 
may transform Eq. (|3^) to a differential equation of third order, after which one may 
perform the averaging over the random magnetic field by assuming isotropy and using 



the spectral density representation for the magnetic field correlators. One finds |^ 
that the probability for helicity change is given by 

Py^^y^[t) = — ^[1 - exp(-z/_Lt)(cosu;t H smcut)] 



Note that in the presence of elastic collisions the spin-flip probability increases and 
tends towards the asymptotic value 1/2 as t — > cx). In the present section we assumed 
that all the collisions take place within a homogenous magnetic domain. Hence, when 
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evaluating the production rate of right-helicity neutrinos, one should calculate the 
probability when t ~ Lq. During this time the neutrino has been subject to a large 
number of independent collisions, each of which have served to adjust the spin-content 
of this state. The scale of the magnetic field felt by the neutrino at each collision is 
given by the free path length L^. Thus we find 

r..^ = ^M^^. (36) 

Assuming z/^ ^ T^r it is then straightforward to deduce a constraint on the product 

fiMTqcD, Lw)<3.5 X 10^ (^^'^ . (37) 

Note that by taking into account elastic collisions one obtains a limit which is more 
stringent than in the coUisionless case. This is due to the fact that spin rotation turns 
the longitudinal part of the spin into transversal, and at each collision the transversal 
part is, in effect, wiped out. This results in a shrinkage of the spin vector and as 
t — s> oo, Pj,^^j,^(t) 1/2. Qualitatively one can also see this in the following manner. 
Let us write the RKE in the familiar form 

^ = V X S - z/^Sx, (38) 
at 

where now V = Vfiz + 2fi^{Bxnx + Byfiy) and the transversal spin is given by S_[_ ~ 
S{y_}^/Vz) with Vj_ = 2fii^B±. It then follows that 

^ = -2u±Sl (39) 

yielding a shrinkage rate S^'^dS'^/dt = — z/_l x 8{H'j_) /V"^, in agreement with Eq. (35). 

The usefulness of the limit Eq. ( PT]) in restricting the neutrino magnetic moment 
depends of course on the magnitude of the primordial magnetic field. Adopting the 
maximum value allowed by nucleosynthesis, given in Eq. (^), and taking Lq ^ Ih 
we see that the tightest possible constraint which in principle can be obtained in this 
manner is /i^<3 x lO'^^fi^- Cosmological and astrophysical constraints on 
fii,, based on the direct production of wrong-helicity neutrinos in photon mediated 
collisions are typically less severe by several orders of magnitude. Thus the presence 
of a primordial magnetic field is a potential bonus for neutrino physics. 

The derivation of Eq. (37) assumes that Lq ^ ^w- This means that a scattered 
neutrino meets always a transversal part of a randomly orientated magnetic field. It 
seems not very likely, however, that the coherence length of the magnetic field at 
T = Tqcd could be as large as the horizon size, especially if the origin of the field is at 
earlier times when the size of the horizon was much smaller. Thus, in the next section 
we consider the corresponding limit in the case of a small scale magnetic field. 
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4 Neutrino propagation in medium with small- 
scale random magnetic field 

An important technical point in the derivation of Eqs. (|3^) and (^) is the procedure 
by which one averages the spin equation of motion over the random magnetic field 
distribution. For the limits (|32D and (0) the exact differential equation derived from 
Eq. was averaged after it was first transformed to a more suitable form. Such a 
procedure is always valid for a regular magnetic field, but in the case of random fields 
we should be more careful. In fact, as we shall now show, in the case of small domains 
with Lq <^ lose there appears aperiodic neutrino spin motion. 
Let us rewrite Eq. ( ^l]) as 

= -4/o*exp[-/,* u^{h)dh] X [Re {H.{t)H^ih)) cos(//^ V{h)dt2) 

(40) 



+Im [H^{t)H^{h)) sin(//^ V{h)dt2)]S,{h)dh • (41) 



Assuming that the collision frequency z/^ depends only weakly on the magnetic field, 
and taking into account that in the leading approximation the potential V{t) as given 
in Eq. (|l]) depends on Bz{t) only while H±{t) are proportional to the transversal 
components, we can average these factors in integrand independently because of the 
isotropy of the system. For neutrinos crossing many small-scale domains with t ^ Lq, 
the size Lq corresponds to a narrow resonance for uncorrelated random fields, as we 
now show. 

To this end, consider the kernel 

K{t~t,) = {HAt)H_{t,))/{Hl) (42) 

with 

(Hi) = {2/3)filB'. (43) 
If the fields are uncorrelatedQ with {B{t)B{ti)) = B^Lo6{t - h) one finds 

^ - hm = ^5{t) . (44) 

For such uncorrelated fields the averaging over of the transversal components then 
results in 6{t — ti) under the integral in Eq. (0) which wipes out the exponent in the 

^Here we consider one-dimensional correlators, but the same result could be obtained by using the 
full 3-d correlator given in Eq. (||). There the factor A is related to the horizon scale because (B^) 
involves integration over all space. 
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integrand and leads (for the initial condition 5*^(0) = —1) to a new damping solution 

S,{t) = -exp(-n) , 

^ il^^^|izE^)).n/2, (45) 

where the damping parameter F is given by 

r = ^fxlB'Lo . (46) 

One should bear in mind that the averaged field B depends on the horizon scale L = In 
and the domain size Lq for uncorrelated fields (see Eq. (^). 

As in the previous section, the spin-fiip probability grows with time. Now we 
assumed that Lw ^ Lq so that we should calculate the probability at largest possible 
time. Setting t = we find 

= ^filB'LoH-'T'^\ (47) 

Note that Tl-*r ^ T^, so that again we should require Tl^r < H at T = Tqcd- This 
results in the bound 

/i,5(TQCD,/H)<6.7x lO-^/i^G (^-^j . (48) 

Note that this bound does not agree with Eq. (^) in the limit Lq Lw- This is 
due to the qualitatively different averaging procedures at large and small Lq. When 
Lq ^ Lw, neither method is reliable. 

Substituting Eq. (|^) to Eq. ( ^HD we can rewrite the BBN constraint on the Dirac 
neutrino magnetic moment as 



1.7 X 10-21/i^ 



^^^ ,.L..^^^.,.S,L.^ ^ (49) 



where the minimum domain scale L'^^^{Tqcd), the left neutrino total collision rate 
Tw{Tqcd) as given by Eq. (pi|), and the Hubble parameter if (Tqcd) are the functions 
of the temperature Tqcd = 200 MeV x T2oo- 

As we have discussed in section 2 there are two scenarios for the choice of the min- 
imum scale L™"(T). Let us first assume that the relic field is the seed for the galactic 
magnetic field, so that the domains with a size Lq > L^^^{T) = 10^ cm(MeV/T) 
survive after the recombination time. In the second scenario^ , with the use of the 

^The first scenario is not necessary since there are other possibihties for the seed field creation in 
the MHD-dynamo theory of the galactic magnetic fields. 
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relativistic plasma dissipation length (|T6]), the domains with sizes Lq > L™" = 
10~^ cm(MeV/T) survive at any relativistic temperature T>2 MeV, but dissipate 
even before the BBN temperature T ~ 0.1 MeV. This last fact does not matter for 
the neutrino spin- flip that populates wrong helicity states mainly around T ~ Tqcd- 
Thus, we find from Eq. (49) two upper limits for the Dirac neutrino magnetic moment: 

1 n-22+6p,, 

for the first case with the relic seeding of galactic magnetic fields, and 

^ 3.2xl0-^°+^^^/i^ 

f^'^^ ' (51) 



{5A)PTi 



•p+2 
200 



for the second scenario with dissipation of the random fields before the BBN time. 
As an example, for the index p = 1/2 we obtain from Eqs. (|50D and (|5TD very re- 



strictive constraints on the Dirac neutrino magnetic moments, /ij,<4.3 x 10 HB/Tr 



20,,_ /T^5/2 



200 ) 



or fi^<4:A X 10~^^fiB/T2oo, respectively. These numbers are deduced by requiring 
that there should not occur full equilibration of one right-handed neutrino species. 
If we require that AN,y<0.1, we should multiply the upper limits above by a factor 
(AN^)^^^"^ ~ 3. However, we wish to emphasize that the very stringent constraints 
on the Dirac magnetic moment above are very sensitive to the model of the primor- 
dial magnetic field. We should also point out that the both scenarios are based on 
the common assumption that the magnetic coherence length is much larger than the 
interparticle distance, Lq ^ T^^, an assumption which is natural considering the 
macroscopic nature of the "glueing" of the magnetic field force lines on hot plasma. 

The remaining issue is the validity of the RKE Eq. (0) when also inelastic colli- 
sions are taken into account. This work is now in progress. 



5 Discussion and conclusions 

We have found out that the general constraints on the primordial magnetic fields, as 
implied by nucleosynthesis, do not exclude the possibility of very tight limits on the 
Dirac neutrino magnetic moments in the presence of magnetic fields. These can be de- 
rived from the requirement that the right-handed components should not be in thermal 
equilibrium at time of nucleosynthesis. We found, however, that nucleosynthesis and 
the dynamo origin of the observed galactic magnetic fields are not compatible with a 
magnetic field model consisting of uncorrelated cells with a scaling index p = 3/2, and 
only marginally compatible if p = 1. If the field reversals observed in the Milky Way 
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are due to the dynamo, then only the scahng law p<l/2 is appropriate (this naturally 
includes the constant background p = as in the case of the Savvidy vacuum ||^). 

The particle physics aspect of the magnetic moment constraints is straightforward: 
one only needs to know the neutrino collision rates just above the QCD phase transi- 
tion. Neutrino spin evolution is then determined by a general RKE. Here we included 
only the effects due to the elastic scattering. Generally speaking, we should also ac- 
count for the dependence of the weak rates on the strong magnetic field. For the spin 
collision integrals this problem is now in progress, together with the generalization 
of the RKE (|34[) to the case of inelastic collisions. It is not obvious how ineleastic 
collisions affect the evolution of the neutrino spin. 

The actual limits on the neutrino magnetic moments depend on two unknowns: 
the strength of the field aX T k, Tqcd? and the size of the homogenous magnetic field 
domain Lq. The magnitude of the field depends on the mechanism by which it was 
first produced, at the electroweak phase transition or earlier, and on its scaling when 
averaged over several domains, if Lq <^ ^v/- ^ given model of the primordial 
magnetic field both these can be estimated, but here we adopted a phenomenological 
view and simply assumed B ~ (Lq/L)^. The magnitude of Lq is a more complicated 
issue as it involves magnetohydrodynamics in the hot plasma of the early universe. 
It seems obvious, though, that Lq ^ 1/T as the build-up of the magnetized plasma 
requires a large number of particles, but that Lq is (much) less that the horizon 
size. One possibility would be that the collision frequency of charged particles (i.e 
conductivity) plays a decisive role in forming homogenous regions. This remains an 
open problem. In the present work we studied the cases with Lq <^ and Lq ^ 
which lead to qualitatively different evolution equations. We also gave some examples 
of the possible order of magnitude of the upper limits on /ij, which could turn out to 
be as restrictive as 10"^°/!^. The limit also depends on whether the primordial field 
provides the seed for galactic dynamo, and thus in the absence of a more detailed 
knowledge of the dynamics of the primordial magnetic fields no definite statement 
about the actual constraint on neutrino magnetic moments can be made. We may 
however conclude that potentially the primordial magnetic field constraint on Dirac 
neutrinos could be very important. 
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